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Abstract 



. We find a static solution to Einstein's field equations on a five-dimensional 

orbifold with a compact S1/Z2 fifth direction and Poincare invariant 3 + 1 
sections. The solution describes a theory with bulk cosmological constant and 

I 3-branes at the orbifold fixed points which carry matter density and pressure 

CN ■ in addition to tension. The radius of the fifth dimension is determined by the 



matter content of the branes. The ratio of the space and time components of 
the metric depends on the fifth coordinate. Thus, the speed of propagation 
of massless fields is path dependent. For example, bulk and brane fields 
propagate with different speeds. 



I. INTRODUCTION 



A solution to the hierarchy problem has been proposed in which the observable universe 
is a 3-brane at an orbifold fixed point of the non-factorizable geometry 

ds 2 = e-^ri^dx^dx" - dy 2 . (1.1) 

The orbifold has fixed points at y — and y = y c where there are branes with tension 
Vo = 12k and V c = —12k, respectively. We will refer to this as the RS model jlj. However, 
the dynamics does not determine the value of y c , leaving it a free parameter. A solution to 
this so called "radion stabilization problem" has been found by adding a bulk scalar field, 
that is, one that has five- dimensional dynamics, to the model ||. Alternative stabilization 
mechanisms have been proposed in 0. 

In this paper we present a static solution of the field equations when the branes carry 
matter density and pressure in addition to tension. The solution fixes the radius of the 
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fifth dimension. To see how this works let's first recall why there is a radion stabilization 
problem in the RS model. The warp factor in Eq. (|1 . 1|) , A = —k\y\, has a constant slope and, 
therefore, the jump of the slope across the fixed point is —2k regardless of where the jump 
occurs. Now if the warp factor had a non-constant slope A'(y), then the jump — 2A'(y c ) does 
depend on y c . Therefore the radius of the fifth dimension has to be chosen to accommodate 
the jump in energy density on the brane. 

Adding matter to the basic setup of RS is necessary for a description of cosmology of 
the model. The cosmology of brane models has been investigated in a number of papers. 
A general formulation was given in Ref. ||. The work in Refs. [§|-[7| is concerned with 
the cosmology of brane models of the Randall-Sundrum type ][]]. These papers recognize 
that the jump in the warp factors across the branes as implied by Einstein's equations is 
enough to give the evolution of the scale factor in a FRW description of the cosmology on 
one brane. Therefore there has been little interest in the behavior of the metric in the bulk, 
see however Ref. [0. It turns out that investigating bulk solutions can yield interesting 
results. The solutions we find exist only for a particular value of brane separation, which 
is determined by the matter content of the branes. Depending on the kind of matter such 
solutions can be either stable or unstable under small perturbations. The salient feature 
of our bulk solutions induced by brane matter, similar to the solutions in Ref. H], is that 
the space and time components of the metric are not identical as they are in the RS model 
without matter. This means that the speed of propagation of massless fields depends on 
their trajectory. For example, bulk gravitons can propagate at a speed different than bulk 
photons. Finding out the strength of gravity and masses of bulk and brane fields in such 
complicated backgrounds is straightforward using the results of Ref. ||. 

In Refs. [|-|^] it is found that the scale factor a satisfies 

d\ 2 /a\ 1 , 1 



where V, p and p stand for the tension, matter density and pressure on the brane respectively. 
We have set the five-dimensional gravitational constant to unity. Our static solution is 
not in contradiction with this brane equation because the solution demands that the three 
parameters V, p and p are such that the right hand side vanishes. 

We present our solution in section [TT]. We study non-static solutions as small perturba- 



tions to our solution in Sec. Ill and present our conclusions in section IV. 



II. A STATIC SOLUTION WITH MATTER 

We denote the coordinates of spacetime by x A , A = 0, . . . , 4, and often use t = x° and 
y = x A . The fixed points are at y — and y = y c . The class of spherically symmetric metrics 
we study is parameterized by three functions of t and y only 0] 

ds 2 = G AB dx A dx B = n 2 (t, y)dt 2 - a 2 (t, y)dx 2 - b 2 {t, y)dy 2 . (2.1) 

Fixing y = (y = y c ) we see that the metric gives a flat FRW cosmology on the brane with 
scale factor R (t') = a(t(t'),0) (R c (t') = a(t(t'),y c )) where dt' = n(t,0)dt (df = n(t,y c )dt). 
We will denote by g^, with p,u = 0, . . . ,3, the induced metric on the brane. 
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The action is 

S = J d 5 x VG [-R - A] + J d 4 x y/=g [-V ] y=0 + J d 4 x yf=g [-V c 



y=yc 



(2.2) 



The constants A, Vq and V c represent the cosmological constant in the bulk (5- dimensional 
space) and on the branes at y = and y = y c , respectively. In addition there is matter 
density and pressure on the branes, introduced into the field equations directly by their 
contribution to the energy momentum tensor: 



T 



AB 



cAB oAB 

f AB + ^ r 5(y) + ^ r S(y-y c ), 



(2.3) 



where T AB is derived as usual by varying the action with respect to the metric, and S AB are 
contributions from perfect fluids of density p and p c and pressures p and p c on the branes, 



Sb = diag (p, -p, -p, -p, 0). 



Einstein's equations are 



R AB _ l _ G AB R 

2 



AB 



(2.4) 



(2.5) 



Here R AB and R are the Ricci tensor and scalar. The gravitational constant is k and we 
work in units of K — 1. 

For the particular metric ([2.1|) Einstein's equations are 
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(2.6) 

(2.7) 
n'b'\ 



nb J 
(2.8) 

(2.9) 



Here a dot is a shorthand for d/dt and a prime for d/dy. The first four equations correspond 
to the 00, 04, 11 and 44 components of Einstein's equations. Conservation of the stress- 
energy tensor gives 



p + 3-(p + p) = 
a 



(2.10) 



on each brane. 
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We look for static solutions to these equations. One may then reparametrize the fifth 
coordinate to enforce 6=1. The resulting equations are most easily solved in terms of the 
warp factors A and N defined by 

a = exp(A) n = exp(iV). (2.11) 



For convenience we also introduce k = y — A/12. In the bulk Einstein equations reduce to 

A" + 2A' 2 = 2k 2 , (2.12) 

2 A" + 3A' 2 + 2A'N' + N" + N' 2 = Qk 2 , (2.13) 

A' 2 + A'N' = 2k 2 . (2.14) 

The solution is straightforward. The metric in the bulk is given, for y > 0, by 

a 2 = alcosh[2k{y-y*)], (2.15) 

n 2 = nl sinh[2% - y*)] tanh[2% - yj\, (2.16) 

b 2 = 1, (2.17) 



where y*, and a* are constants. For related bulk solutions see Ref. [10]. 

To complete the solution we must examine the field equations at the brane. The delta 
functions on the right hand side of Eqs. ( ^.6[ ) and ( |2.8j ) must be saturated on the right hand 
side by the second derivatives. To this effect we demand 

(2.18) 
(2.19) 
(2.20) 
(2.21) 

The solution to these are two equations of fine tuning, 

144A; 2 = (Vo + 2p )(V - Po - 3 Po ) = (V c + 2p c )(V c - p c - 3p c ), (2.22) 

and two equations fixing the parameters y* and y c : 

I2kta,nh[2ky # ) = V + 2p , 12ktanh[2k(y c - y m )] =V C + 2p c . (2.23) 

The first two conditions, Eq. ( |2.22 ) are similar in spirit to the two fine tuning equations 
in the RS model that set the square of the tension in each brane in a fixed proportion to 
the bulk cosmological constant. By comparison, the stabilization mechanism of Refs. p|JTT 



2 A \y=0+ 


= —=Vq - Ipo, 
o 3 


—9 A'\ 

AS\ \y=y c - 




o 3 


2(2A' + N')\ y=0+ 


= --Vq+Pq, 


2(2A' + N')\ y=yc _ 


= ~^V c +p c . 



requires only one fine tuning, which is equivalent to setting the cosmological constant to 
zero. The last two equations describe what was alluded to in the introduction, that in order 
for the warp factor to jump by the appropriate amount the location of the brane must be 
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chosen accordingly. Hence the radius is stabilized, or, should we say, equilibrated (stability 
is investigated in the next section). 

Note that positive tension branes (V > 0) can only exist for p < —p. If ~p > p > — ~p 
on both branes, solutions exist only for V < —12k and V c < —12k. But if p > |p, V can 
be larger or smaller than —12k. 

If we insist, as in the RS model, that Vq = V 2 = 144fc 2 then the conditions (|2.22 ) are 
identical to the static solution conditions in Eq. ( |1.2| ). We recover the RS model taking 
vanishing matter density and pressure and the limit — > oo. 

However we need not insist on imposing Vq = V 2 = 144fc 2 . In fact, we are now free to 
chose any tension provided the conditions fl2.22|) are satisfied and a solution to Eqs. ( |2.23| ) 
with y c > can be found. 



III. SMALL PERTURBATIONS 



Armed with the new solutions with static matter density, we proceed to investigate the 
time dependence of small matter perturbations. Let us denote the static solution of the 
previous section by n = e$ , a = , and b = 1. We look for solutions to the field 
equations, Eqs. (|2l 



(3.1) 



M, of the form 

n = n (l + 5n), 
a = ao(l + 5a), 
b = b (l + 5b). 

In addition we set the density on the branes to po + <5po and p c + 5p c and the pressure to 
p + 5po and p c + Sp c . 

We count orders of the perturbative expansion parametrically in 5p and 5p. That is, we 
re-scale 5p —> e5p, count powers of e and set e = 1 at the end of the calculation. In particular 
this implies that we make no assumption as to the relative importance of temporal or spatial 
derivatives [[]]. 

To derive the linearized equations in the bulk, we use the parameterization in Eqs. ( p.l| ) . 
The 00, 04, 11 and 44 components of Einstein's equations give 

5a" + A' {45a' - 5b') 

d 



-((N^-A' )5a + A' 5b-5a') 
25a" + 5n" + (QA' + 2N' )5a' + 2(A' + N' )5n' 

-{2A' + N' Q )5b' -\ (25a + 5b 

Tin 



(2A' + N^)5a' + A' 5n - —5a = 
The solution to these equations gives 5b and 5n! in terms of 5a: 



Ak 2 5b, 
0. 



(3.2) 
(3.3) 



12k 2 5b, 
Ak 2 5b. 



5b 
5n 



A, 
l 

An 



n. 



;5d - (2A' + N^Sa' + Ak 2 5b 



(3.4) 
(3.5) 

(3.6) 
(3.7) 
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In Eq. ( p.7|) 5b is understood as shorthand for the solution of Eq. ( |3.6| ). In Eq. flOf ) F is a 
function of y satisfying 



F' — 



A% - Ak 2 

An ' 



0. 



(3.8) 



It must be observed that F may be discontinuous at y — 0. In fact, continuity of 5b at y = 
requires F(0+)+F(0-) = 0. The solution to Eq. is F(y) = F*/sinh(4% -y*))> with 
F* a constant. 

We connect the bulk solutions for y > and y < demanding continuity of the fields at 
the brane, y = 0, and using the jump equations for the discontinuous derivatives at y = 0. 
The latter give jump conditions for the perturbations 



o+ 



5a' 

{25a 1 + 5n')\ 0+ 



{^V + p )5b + 5p 
{-\v +p )5b + 5p Q 



Similarly the jump equations at the second brane are 

.1 



- 5a'\ Vc - = - 
- (25a' + 5n% c - = : 
In addition, conservation of energy gives, on the brane 



V c + p c )5b + 5p c 



2 

-W c + p c )5b + 5 Pc 



5p + 3(p +Po)5a\ 



and 5p c + 3(p c + p c )5a\ 



y=yc 



0. 



(3.9) 
(3.10) 

(3.11) 
(3.12) 

(3.13) 



The right hand side of these equations could be a non-zero constant. However, we set it 
to zero since we are not interested in constant shifts in the mass density (since these are 
accounted for in the exact solution of Sec. ||). 

The jump equations at y — 0, Eqs. ( jO|) -( [37L0l ), determine 5a' in terms of 5a and give 
an equation for 5a, namely 



1 A'r 



j6d + 2~W 5P ° ~ « 5p ' 
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(3.14) 



The time dependence of the matter is fixed by the continuity equation, Eq. ( |3.13|) . Given 
an equation of state for the perturbations, 5po/5po = Wo, one can solve this equation. Let's 
rewrite the equation as 



5a — Tl5a = 0. 



(3.15) 



Then the coefficient Tq can be expressed in terms of the initial parameters. In terms of 
t = (Vb + 2p )/12k = tanh(2%*) we find 



2k 2 n\y = 0) 



- (3u + l)t c 



(3.16) 



Hence, Tq > for equations of state with wo < 1/3. But for u> > 1/3 rjj can be either 
positive or negative depending on the value of t with u = |(4 — 1) being the dividing line. 

We see that the time dependence of 5a on the y = brane is exponential for oo < 1/3, 
but can be oscillatory for u > 1/3: 

Sa\ y=0 = c e r °* + d e- Vo \ (3.17) 

where T is either real or purely imaginary. Substituting this solution in the equation for 
8a' gives 

5a'\ y=0 = 0. (3.18) 

This implies that the jump in F vanishes. Therefore F(y) = 0. 

An entirely analogous solution is obtained on the brane at y — y c by replacing the 
subscript "0" in Eqs. ( fj. 15| ) ( P- 17[ ) by the second brane subscript "c." 

To extend the solution on the branes into the bulk we must fix the gauge. There is a 
gauge freedom, that is, reparametrization invariance consistent with the form of our metric. 
Starting from the metric 

ds 2 = n 2 (t\y')dt' 2 - a 2 (t',y')dx 2 - b 2 (t' \y')dy' 2 } (3.19) 

we look for infinitesimal transformations 



t' = t + T(t,y) (3.20) 
y' = y + Y(t,y) (3.21) 

that leave the form of the metric invariant and has fixed points at at y = and y = y c . Here 
T and Y are infinitesimal. The only constraints on these functions come from the absence 
of off-diagonal terms in the metric, 

n 2 T' - b 2 Y = 0, (3.22) 

and from the fixed points, 

Y(t,0)=Y(t,y c ) = 0. (3.23) 

Under the gauge transformation the metric variations are 

A5n = NqY + T, (3.24) 
Ada = A' Y, (3.25) 
A5b = Y'. (3.26) 

For simplicity we have indicated the variation about a static solution with bo = 1 and 
do = ho = 0. We would like to use this gauge freedom to impose the gauge condition 

5b(y,t) = 0. (3.27) 
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However, this cannot be done in general since there is only one constant of integration in 
the gauge condition ( p. 26 ). Instead consider covering the space with two different patches. 
We label the metric perturbations 8a\, 8n\ and 8bi on the patch defined by < y < y c /2 + e 
and by 8a 2 , 8n 2 and 8b 2 on the patch defined by y c /2 — e < y < y c . Given a solution with 
metric 8a, 8n and 8b we choose the new coordinates by choosing 

Yy = - F 8b{y,t)dy and Y 2 = - F 8b(y,t)dy. (3.28) 

JO Jy c 

Thus 8b\ = for < y < y c /2 + e and 8b 2 = for y > y c /2 — e. 

The solution for 8a\ and 8a 2 follows immediately from Eq. (|3.6|) setting 8b = 0: 



5a x (y,t) = F #1 /4k + A' ^(t) and 8a 2 (y,t) = F* 2 /4k + A'^ 2 (t), (3.29) 

where are arbitrary functions of t but independent of y. We can now use our jump 
conditions to specify these completely: 

A' (y) A' (y) 

5ai{y,t) = A' °(0+) 6a ^ =0+ 5a ^y^) = Jih _} Sa \v=vc- ( 3 - 30 ) 

In the overlap region, \y—y c /2\ < e, these solutions are related by a gauge transformation: 

Sa 2 - 8a x = A' F" Sb(y, t)dy. (3.31) 
Jo 

Thus we obtain 



«<= „ 8a 



6b(y,t)dy 
o A' y=yc- A' Q 



8a 



y=0+ 



c c e Tct + 4e~ rci + c e Tot + d e- rot ^ 



kta.nh(2k(y c — y*)) fctanh(2%*) 
where it is understood that r 0iC can be purely imaginary if u 0jC > | 



IV. DISCUSSION AND CONCLUSIONS 

In Sec. m we gave a class of solutions to the field equations with matter on the two 
branes. The solutions are static. The price to pay for time independence is a fine tuning of 
the amount of energy on each brane, as expressed in Eq. (|2.22|) . These fine tunings are no 
worse than the corresponding fine tunings in the RS model. 

The salient feature of the solutions found in Sec. |I| is that the physical size of the fifth 
dimension is fixed. This suggested the exciting possibility that the radion stabilization 
problem is not an issue in models with matter on the branes. 

However, static solutions are not cosmologically acceptable. Nevertheless, one wonders 
if even for non-static solutions with matter the radius is stabilized. Short of finding an exact 



non-static solution we have exhibited in Sec. |TJ an approximate solution by linearizing the 
field equations around the static solutions. 
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The result of the linearized analysis is summarized by Eqs. (|3.17Q and ( |3.32j ). The first one 
establishes that the metric perturbations have exponential time dependence for equations 
of state with uj — p/p < 1/3, but can be oscillatory for uj > 1/3. The second gives the time 
dependence of the radius of the space. Indeed, a measure of the radius is 



^-G M Ndx M dx N (4.1) 



along a line of constant t and x: 



L = y c + I*" 8b(y,t)dy. (4.2) 
J o 

For it) < 1/3 the radius grows exponentially, at least while the exponential is small enough 
that the linearized solution remains a good approximation. On the other hand, for uj > 1/3 
the radius can oscillate about the equilibrium value y c or can grow exponentially depending 
on the value of t = ^±f^. 

One gains some understanding of the behavior of the scale factor a on the brane by 
considering the full, non-perturbative description of its evolution at an orbifold fixed point. 
Since we do not wish to insist that V 2 = 144/c 2 we modify Eq. ( |1.2| ) to allow for an uncon- 
strained tension, 

+ (!)^^- 3 ")-^ +3rt+ ^ 2 - 2p - < 4 - 3) 

Let the equation of state be p = wp. Using the equation of conservation, Eq. ( |2.10|) , one 
has p = poa^ 1+w ^ / a 3 ( 1+w \ where po and a are the density and scale factor at a fixed time. 
Then one can rewrite the equation for the form factor as the equation for a particle with 
displacement r = a 2 in a potential, r = —U'(r), with 

U{r) = dr 2 -'^ 1 ^ + C 2 r 2 - 3 ^ + C 3 r 2 , (4.4) 

where 

3(1+«0t/ 

Cx = -**2 L 

18 V ' 

„2 6(1+™) 

C, = -^g— (4-6) 

C 3 = 2k 2 -^ (4.7) 

Recall that for u > — 1 our solutions must have negative tension branes V < 0, so C\ > and 
C 2 < 0. For !>u;>— |, V < —12k and, therefore, C 3 is negative and the only extremum 
of U(r) is a maximum. But for w > |, V can be larger or smaller than — 12k and, therefore, 
the sign of C3 can be either positive or negative. The extremum of U(r) can, therefore, be 
a maximum or minimum, in accordance with our analysis of small perturbations. 

The physical interpretation of the solutions discussed here is not straightforward. One 
cannot, as in the case of the RS model, simply renormalize fields on either brane to cast 
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their kinetic energy term in standard form, since n ^ a. Moreover, a null trajectory parallel 
to the branes has dx/dt = n/a. Since n/a = (n*/a*) tanh[2/c(y — y*)] differs between the two 
fixed points, one can arrange for superluminal signal travel. For example, on one brane one 
can send a graviton across to the other brane, relay the signal along the brane by photons, 
and then relay the signal back to the first brane via gravitons. Neglecting the time of travel 
between branes, we see that the interval between emission and reception can be shorter than 
the time for a photon to travel between the same two points directly on the first brane \Y1 . 
A complete discussion is the subject of a forthcoming paper. 

Since for standard fluids in equilibrium the equation of state must have uj < 1/3, the 
exponential growth of the scale factor and radius make this model an unlikely candidate for 
cosmology. However, it may be possible to stabilize the radius by modifying the model, say, 
by adding matter in the bulk, e.g., scalar fields 0]. The cosmology of such a model may 
indeed be acceptable ||. 

Note Added While this article was being completed three related works have been sub- 



mitted to the archives. Ref. |L3| discusses general solutions with different space and time 
components of the 5-d metric. We have not attempted to check if an explicit coordinate 
transformation can convert our bulk solution to a form used in Ref. ||13|| . The authors also 



discuss Lorentz symmetry violations due to different propagation speeds of bulk and brane 
fields. Ref. [I4"| uses super symmetric gauge dynamics for stabilizing the radius and Ref. [15 



uses the Casmir force due to a bulk scalar field. 
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